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The cosmological primordial power spectrum is known to be one of the most promising observable 
to probe quantum gravity effects. In this article, we investigate how the tensor power spectrum is 
modified by loop quantum gravity corrections. The two most important quantum terms, holonomy 
and inverse- volume, are explicitly taken into account in a unified framework. The equation of 
f*"^ . propagation of gravitational waves is derived and solved for one set of parameters. 
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I. INTRODUCTION 



The inflationary scenario is currently the favored paradigm to describe the first stages of the evolution of the 
Universe (see, e.g., [H for a recent review). Although still debated, it has received many experimental confirmations, 
£~ ' . including from the WMAP 5-year results Q, and solves most cosmological paradoxes. 

On the other hand, a fully quantum theory of gravity is necessary to investigate situations where general relativity 
(GR) breaks down. The big bang is an example of such a situation where the backward evolution of a classical 
s ! ' space-time comes to an end after a finite amount of time. Among the theories willing to reconcile the Einstein gravity 
^ | with quantum mechanics, loop quantum gravity (LQG) is appealing as it is based on a nonperturbative quantization 
of 3-space geometry (see, e.g., [|[ for an introduction). Loop quantum cosmology (LQC) is a finite, symmetry reduced 
model of LQG suitable for the study of the whole Universe as a physical system (see, e.g., Q). 
. In this article, we consider the influence of LQC corrections to general relativity on the production and propagation 
of gravitational waves during inflation. We first derive the equation of propagation of gravity waves with both 
holonomy and inverse-volume corrections. This equation is then reexpressed with the commonly used cosmological 
■ variables. It is finally solved for a specific set ofparameters and the primordial power spectrum is derived. The aim 
of this work is to conclude our previous studies Q and where, respectively, only holonomy and only inverse- volume 
r* ■ corrections were considered. By combining both terms, we show that the inverse-volume correction dominates over 
the holonomy one and dictates the overall shape of the tensor spectrum. 

Quite a lot of work has already been devoted to gravitational waves in LQC 0. Our approach assumes the 
background to be described by the standard slow-roll inflationary scenario whereas LQC corrections are taken into 
account to compute the propagation of tensor modes. This approach is hcuristically justified (to decouple the physical 
effects) and intrinsically plausible (as, on the one hand, the LQC-drivcn super-inflation can only be used to set the 
proper initial conditions to a standard inflationary stage if the horizon and flatness problems are both to be solved 
Q and as, on the other hand, it seems that the quantum bounce can trigger on a standard inflationary phase @). In 
addition, very few studies so far have taken into account both the holonomy and the inverse-volume corrections. This 
latter term is somehow more speculative than the former one as it was shown to exhibit a fiducial cell dependence 
(see, e.g., [l(|). For the sake of completeness it is however obviously worth considering the fully corrected propagation 
of gravitational waves. 
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II. EQUATION OF PROPAGATION FOR THE GRAVITON 

The derivation of the equation of propagation of gravitational waves with both holonomy and inverse- volume 



corrections extensively uses the material developed in jll| : notations, conventions and framework of this work are the 
same and will not be explicitly restated. We begin by considering a Fricdmann-Lemaitre-Robcrtson- Walker universe 
with a spatial metric q a b which will be perturbed to account for gravitational waves. Hereafter, N and N a are 
respectively the lapse function and the shift function. The metric components read as follows: 

.900 = -N 2 + q ab N a N b = -a 2 ( V ), (1) 
5oa = aWV b = 0, (2) 
9ab = q a b = a 2 {rii){8 ab + h ab ). (3) 

As usual in the formalism of LQC, we use the Ashtekar variables for an homogeneous and isotropic background: 
the connection A l a , and the triad density Ef . They can be written as a function of two other variables (k,p) as 

E? = pSf 

K = ksi, r* = o, 

N a = 0, N = ^p. (4) 

Hamilton- Jacobi equations will be used to determine the perturbed part of the Ashtekar variables. The Hamiltonian 
constraint reads as 

H[N] = hj/* XN \ detE \~ hE i Ebk{eijk ^ (5) 
where F l ah = d a A\ — dbA l a + e n A^A^ is the field strength. The Hamiltonian for a matter field $ is given by 

^2 i rpa rpb 




i P i + E?E?d a $d b $ 



detE^ 



Hmatter = / <Px \ f ' ° + W ' \detEZ\V '(*) . (6) 



With Eq. (j4j) and these Hamiltonians, the background is described by 



2k 
and 

,2 



H f G ond [N] = ±- I d 3 xN [~6Vpk 2 ] , (7) 



H ma tter[N] = I X ( \*±+$V<?) ) . (8) 



Perturbing the canonical variables (and going through the appropriate Poisson bracket) leads to: 

k 2 



H G [N] = ±-J d 3 xN 



2k 1 



(9) 



where only the tensor perturbations (i.e. gravitational waves) are considered in 5Ef. 

This classical Hamiltonian is to be modified by quantum corrections. Because loop quantization is based on 
holonomics, i.e. exponentials of the connection rather than direct connection components, one needs to substitute in 
the gravitational sector 

k -+ Sin(m/27 ^ (10) 

mfi'y 
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where p, is a new parameter related to the action of the fundamental Hamiltonian on a lattice state. In addition, 
because of inverse powers of the densitized triad which, when quantized, becomes an operator with zero in the discrete 
part of its spectrum, the matter and gravitational Hamiltonians must be modified by introducing the function 



( 1 ~y 



(11) 



At a semiclassical level, i.e. g « 1, the same parametric form of a can be used in both the matter Hamiltonian 
and the gravitational Hamiltonian. However, the two positive and real valued constants A and n may differ from one 
sector to another. In the following, (S, A s , s) and (D, A^, d) will therefore denote (a, A, n) for the gravitational 
sector and the matter sector respectively. With these two corrections, the Hamiltonians read 



Hg[N] = i- / d 3 c^ 



-6VP 



sin \v)k 
M7 



n 2 



d 6 x[ -D(q)^+ptV(<5>) , 



(12) 
(13) 



with Hq^ the effective gravitational Hamiltonian describing the homogeneous background. The equations of motion 
for (fc,p), i.e. the background equations, can be obtained in the Hamiltonian formalism 



P = {p, H C J> [N] + H matter [N]} ; k = {k, H e J f [N] + H n 
sin{2p, r yk) \ 



leading to: 



AN]}, 



k = 



2-p- S{p, SE) 

^ 3ffmatter 

Wo dp 



2/2 7 J 



dS f sin(fi'yk) \ 2 S 



dp 



M7 / 



_ 2 — 

/ sin(ft , yk)\ ^ ^ 9 f sin(ft'jk) 
V M7 / dp\M 



(14) 

(15) 
(16) 



The same modification is applied to the perturbed gravitational Hamiltonian. Denoting HQ hen the effective per- 
turbed quantum-corrected gravitational Hamiltonian, it reads with both holonomy and inverse-volume corrections 



rjPhcn 

G 



[N] 



1 

2/7 



d 3 xNS{p, 5E?) 



/ sin ii^k\ 




( sin [i"ik\ 




2p 3 / 2 


\ M7 ) 



P 



3/2 



{8E*8E d k 5%) 



(5 cd P k E^d e d f E%) 



We now turn to the equation of motion of the graviton. The perturbed densitized triad is 

SE9 = -lphl 



(17) 



(18) 



As has been done for the homogeneous canonical variables, it is possible to define the equation of motion for the 
perturbations: 



SEf = {6E?,H£ hen [N} 

= -{5Kl(x),SEUy)} 



HmatterW]} 

6 



5(8K>) 



(H^ n [N] + H matter [N]), 



SK* = {SK,H^ en [N} + H matter [N}} 



{8Ki{x),5E b Ay)} 



S(SE) 



{Hq [N] + H. n 



AN]). 



This leads to: 



SE ? = --(ph«+ph1) 



= S(p,5E) 



p-SKl.8l.5\-{^^).5Et 



(19) 
(20) 
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By combining those equations and using the expression of p, one obtains the expression of 8K l a as a function of h l a 
and of h\. The expression of 6K % a is: 



6K « = 28* ' 2 



1 / sin(2^i'jk) \ 



2/2 7 J »■ 



(21) 



The equation of motion will lead to another derivative with respect to r\. The Hamilton-Jacobi equation for the 
perturbed connection can now be used to find the final equation of propagation for gravitational waves: 



SK* 



S S 2 d V a+ \ 2/2 7 / ° ° ' d V \ 2/2 7 J 



{SKi(x),SE b Jy)} 



S(SE b ) 



(H£ hen [N] + H„ 



As 



§ H Phen 



2 k 



d 3 {x) ■ N ■ 



5S 



S(SE) 



2 k 



d 3 (x)NS 



2 

2pl 



2 f S in(2fj 7 ~k) \ ■ 2 



-{8 bd -V k -8 eS d e d f {5E d k )) 



where [...] stands for the term beginning with[— 
V 2 (SE$) = -±p-V*hi) 



sin f-t^k 
£7 



.] in (jnj, one obtains (with 5 e ^d e d f (SE 



J <5(<5^) 



-5 



1 ^ sinQjrfk)^ 2 hi _ ^ sin(2fxjk) ^ [K + ( sin(2^k) ^ ^ 



2/27 / \ 5 



2/i7 / *, 



V 2 M 



s s 2 dr] a + \ 2/27 ; a a ^v 2/27 ; 



.[JV] 



(22) 
(23) 



After quite a lot of algebra, the equation of motion of the graviton can be derived: 



2S 



sin(2fj,jk) \ ^ ( 
2^7 / ° \ 



p_dS_ 
S~dp 



S 2 V 2 K + S'ToK 



(24) 



where 



^dpj V M7 / 



IT 



1 dH matter ( SE^Jl 



Wo dp 



cos{2fi r yk) 



5H n 



6(6 E?) 



A' 



1 5S 
2 Vp 6(6Ef) 



OS / sin(fi'jk) 



p-—cos(2fi'jk) 
op 



/O 



hi. 



As usual, requiring an anomaly-free constraint algebra in the presence of quantum corrections requires A l a to vanish. 
It should be noticed that the inverse- volume correction is involved in each term, through the S and D factors, whereas 
the holonomy correction is only involved in the h\ term, in Tq and in IIq^ . 
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It is worth studying a bit more into the details of this Hq source term as it seems to have been misunderstood in 
several works. In particular, it has often been either neglected or miscomputed. Without holonomy and inverse- volume 
correction, this term reads as 



IT 



1 dH mat 
3Vb dp 



SH mat 



6(5 Ef) 



with, in this case, 



At the zeroth order in gravitational perturbation, one can show that 

r,2 



Elmatte 

and the nonlinear H matter is given by 

H-matter — Hmatter 

thus leading to 

5H maU , 



d xN 



1_ 

2 pi 



+p*v{4>) , 



(25) 



(26) 



(27) 



2p 3 



nse 



3 ^ i\$~ v ^ 



(28) 



(29) 



6{SEf) 2 y/p 

Restricting to the first order in perturbation, the derivative with respect to p can be evaluated and one finally obtains 

6EJ6l6l_ SH n 



3V dp 



P 



S(SEf) 



(30) 



This easily establishes that classically II„ = 0. However, when LQC corrections are taken into account the source 
term may not vanish anymore (because of the derivative of D with respect to p for the inverse- volume correction and 
because of the cosine term for the holonomy one). 

When only inverse- volume corrections are considered, the source term is still given by Eq. (|25|) but the matter 
Hamiltonian now reads 



JT , tt(S) 

"-matter T ^ matter 



(31) 



d s xN 



Z r>2 



■p*v(</>) 



^Vp 



5E?5E 5idl\D(p,SE?) 



2p 3 



V{4>) 



which leads, at the leading order, to 

SH matter 

WE?) 



N 



5E° 



SD 



2V? a " \ 2p3 



2pi 6(5 Ef) 



and 



1 dE matter 5E-M -l fsmSiSt 



N- 



3 D 



3Vb dp p 3 

We finally obtain 

ri,(IV) _ 1 / $Ej5 J a 5l \ dH ma tter 5H mat ter _ v\ 



n 



Qa 



dp 



6(6E?) 2 pi 



2 op 2 p2 



1 ( 6E-8j6l \ dD t 8D 



p dp 5(8Ef) 



(32) 



(33) 



(34) 
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However, because of the anomaly-free condition (see Eq. (27) of [HI), this term is vanishing. This means that, at the 
leading order, ITg^ /y ' = 0. 

Considering now the holonomy correction alone, one can expand the cosine term in IIq^ and show that 

lft (Mo) = -2/2 7 sin 2 (^k) 



T i,(holo) _ _ 0n .2 ( r T\ 1 9H maU er / ggjgg 



3Vb dp \ P J 

Considering simultaneously the two types of corrections and using the explicit expression of the matter Hamiltonian, 
one obtains the full LQC source term: 



as expected from the vanishing inverse-volume source term. 



III. SCHRODINGER EQUATION FOR THE FOURIER MODES 

The energy density and the pressure of the cosmological fluid can be written as 

_ 1 8H ma tter _ 1 SH mat t er [N] 

9 ~ V p-i SN ' P ~ ~ NV 6(y/\fetE\)' 



With the Hamiltonian constraint, one obtains 

2k 



= ^1 d*xS 



5H r , 



SN 



(37) 



which finally leads to 



3 S ( sin(jJ/yk)\' 



. (38) 



Defining TL as the Hubble parameter with respect to the conformal time (H = a 1 da{rf) j 'drj)) , we obtain the quantum 
Fricdmann equations: 



H 2 = S 2 



sin(2/j,'yk) \ ' 



2/2 7 J 



5 2 |^fl-/2V?|p), (39) 



which lead, with p c = 3/(k/i 7 p), to 



M a = a" ( S _ P.). (40 ) 



This equation, which has already been found in [12|, includes all the LQC corrections and shows that the holonomy 
term, leading to the bounce, is the most important one as far as the background in concerned. This conclusion will 
be radically modified for perturbations. 

The equation of motion for the graviton can now be reexpressed in terms of the commonly-used cosmological 
variables. By taking into account Eq. (|3"5]l and jx 2 p = l 2 PLl one obtains 

S 2 Tq = 2 ( (W ( = l-f. (41) 
\fj.dpj V / 3p c 

The multiplicative factor of h\ in Eq. (|24|) can be reexpressed as a function of the Hubble parameter 
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Finally, the source term can be explicitly computed 

n' 



Kp_p 

S p c 2 



<I> 2 



1 D a 

D(q)a 2 { 1 ~ 6Dd 



As in [TT| , we use the effective parametrization 5=1 + X s (q) ' with q = (a/lpL) 2 . The equation of propagation can 
now be written as 



with 



(43) 



M 2 (a) = K^-a 2 I \ P - 



& ( IDa 

D(q)a 2 \ ~ 6D a 



p c V 3 

This can be usefully expressed as an equation for the spatial Fourier transform h k of h 



(44) 



h k + 2 C -ll-l-^)h k + (S 2 k 2 - M 2 {a))h k = 0. 
a \ 2a b I 



(45) 



The variables are changed according to 4> k = h k a/\/~S, leading to a Schrodingcr-like equation 



id dS 3 I S \ IS 

> k+{S 2 k 2 -\- + M 2 (a)-~- + -l-\ --- 



= 0. 



(46) 



IV. POWER SPECTRUM 

The main question to address is to investigate if one correction, either holonomy or inverse- volume, dominates over 
the other as far as the production of gravitational waves during inflation is concerned. The system describing the 
dynamics is 



W 2 = a 2 'ip S-^ 

Pc 



= + 2- ( 1 - ) $ fc + a 2 DVM&), 

a \ 2a D ] 



2;„2 



= 4>k + < s z k 



M 2 {a) 



a S 3 I S 
aS + 4 I S 



IS 
2S 



0k, 



which is unfortunately much too difficult to be analytically solved. We therefore turn to the approach developed in 
0,0]. The background evolution is assumed to be classical {D rj 1) with the scale factor given by the usual slow-roll 



approximation a(rj) = Iq\t]\ 1 £ . In this case, the effective Schrodinger equation 
reads, to first order in X s , as 



drf 2 



E h {n) = s 2 k 2 = 



l l -f) 



k 



E k {r))-V{ V ) Ml) = 0, 



(47) 



V( V ) 



2 + 3e 6 1 (l+4e) 2(1 



2 1 n ™ 

if k p c 



(1-6) 



'PL 



12 1 (l+4e) |7?|s _ 2+e(a+2) 

K p c Iq 



+ s(l + 2e)\r,\ s ^- 2 - -s(s - 1 + e(2s - l))|?7| s ( 1+£ )- 2 



(48) 
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To implement initial conditions, we consider the limit 77 — > —00 where the adiabatic vacuum holds. Of course, if 
higher order terms in X s were to be included, the vacuum would not be the same anymore. However, we have checked 
that the adiabaticity condition would still be fulfilled in the relevant wavenumber range. 

It is possible to solve analytically this equation, at least for one set of parameters: s = 2 and e = 0. It becomes 



d 2 (j>k 
dr] 2 



1 + 2A S 



IpL 



1 



3 1 1 



X, 



IpL 



12 1 1 



1 



0. 



(49) 



By some appropriate changes of variables, this equation can be turned into a Whittaker equation. The solution can 
be expressed with Kummer functions and the Wronskian condition (^kd^X — 4^dri4>k = 16iir/M PL allows one to 
normalize the modes. The field is then given at the end of inflation by 



<Mc) 



2v/27T 



M PL (kV2Z) 

and the resulting primordial tensor power spectrum is 

1G 



* +t *U Q+ n-v,l + 2fi,ic 



P T {k) 



M 2 



2/t 



r(fc-i) 



r(o) 



with 




k 2 - Z 1 - 4 



2,2 
PL 

I 2 



Vwzk 2 



'0 



(50) 
(51) 

(52) 

(53) 
(54) 



where Z = (Ipl/Io) 2 ^s and 7 2 = 3/(Kp c l PL ). The ultraviolet limit of this spectrum can be easily derived and leads 
to 



p«v(k) = 16. 3 (^) (i + g(i- 4e) )r 

with uj = 7 Ip^/lo- On the other hand, the infrared limit is given by 

2 



Pi R (k) = 16tt 3 ( ) (Z(l - 4w))-ifc 3 e ,r 



(55) 



(56) 



Those results show that the k — > +00 limit of the power spectrum is in agreement with the general relativistic 
behavior with the addition of a slight tilt. The ultraviolet spectrum is nearly asymptotically scale invariant. This is 
not surprising as both the holonomy correction (encoded in the k~ 3 W term) and the inverse- volume correction (encoded 
in the (l + §p-(l — 4e)) term), taken individually, lead to this behavior. The infrared limit is more interesting as, 
in this case, the holonomy and inverse-volume corrections lead to very different spectra. The result obtained here 
shows that the power spectrum is exponentially divergent, in exact agreement with the limit obtained with the 
inverse-volume correction alone. This proves that, under the standard inflationary background evolution hypothesis, 
the inverse-volume term strongly dominates over the holonomy one. This is to be contrasted with the background 
evolution in the very remote past where the holonomy term alone leads to the replacement of the singularity by a 
bounce. 



V. CONCLUSION 



This work derives the fully LQC-corrected equation of motion for gravitational waves. This equation is expressed 
in terms of cosmological variables and is explicitly solved for a given set of parameters in a standard inflationary 
background. It is shown that the spectrum remains exponentially infrared divergent, as for a pure inverse-volume 
correction. This reinforces the use of primordial gravitational waves as a strong probe of loop quantum gravity 
effects. The next step is naturally to build a fully consistent model which includes all the corrections for both the 
perturbations and the background . 
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